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I. INTRODUCTION 

First introduced by Zel'dovich [jj, the anapole moment (AM) of a quantum system is a quantity 
that involves both the electromagnetic interaction and parity nonconservation (PNC). Ignored for 
a long time, it is not before studies by Flambaum and Khriplovich Ol that the concept acquires 
practical interest. In their work, these authors especially emphasize that the AM would grow with 
the size of the nucleus, making heavy nuclei natural candidates for an observation. Implying the 
hyperfine structure, a first measurement was performed a few years ago in the 133 Cs nucleus 

The deuteron AM has also received some attention recently 0>0 ; lIQ>0- Though its interest 
is largely academic (an experiment is not feasible in a near future), it offers the advantage of a 
laboratory where methods and ingredients relative to an estimate can be studied in details. These 
studies have been concerned with the pion-exchange component of the PNC nucleon-nucleon (NN) 

n q 

force. Calculations were based on assuming an effective-field-theory description [J, la], zero-range 
NN strong forces 0], or more realistic NN strong forces 0- The use of an alternative field-theory 
description was also proposed 0l- The next step concerns the extension of these results to include 
the component of the PNC force due to vector-meson exchange that could contribute as much as, 
if not more than the pion exchange. 

Determining the AM supposes to calculate the effective current that couples to the photon. 
As it involves parity nonconservation, the current has necessarily an axial character, making the 
requirement of current conservation non-trivial. Individual contributions are proportional to the 
weak coupling, Gf, while fulfilling the above property implies that the effective current contains 
the factor Gpg 2 , which vanishes in the limit of a zero momentum transfer. Getting this result 
demands particular care with contributions ensuring gauge invariance. When dealing with vector 
mesons, this task becomes essential. In particular, it has to be done consistently with the PNC 
interaction model that is employed in calculations. Some contributions were given in Ref. 0- I n 
the present work, we intend to complete this work with the double aim to satisfy gauge invariance 
and consistency with the PNC interaction model; DDH potential, given by Desplanques, Donoghue, 
and Holstein, is our case [lfl| . Some estimates of the vector-meson-exchange contributions to the 
anapole moment will be presented in the deuteron case. 

The plan of the paper is as follows. In Section[nJ we first present the various ingredients pertinent 
to the interaction: parity-conserving (PC), parity-nonconserving and electromagnetic (EM) ones. 
We subsequently provide the expressions for the PNC two-body currents at the lowest l/m N order 
and show how they allow one to fulfill current conservation. Section ITTT1 is devoted to applications 
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involving the deuteron. This includes the deuteron description, especially the determination of the 
PNC components; the expression of the anapole matrix elements from both the one- and two-body 
currents; and a numerical estimate in terms of the PNC meson-nucleon coupling constants. A 
discussion of the results is given in Section IIVI This is completed by Appendix El that contains 
expressions of the two-body currents in configuration space. 

II. PNC NN INTERACTION, CURRENTS, AND CURRENT CONSERVATION 

The anapole moment is a special electromagnetic property of a system in which parity conser- 
vation is violated; therefore, our first step in calculation is formulating the EM current operators. 
Throughout this work, we assume for a nuclear system the validity of non-relativistic (NR) limit 
and keep only terms of leading order. 

The one-body term (p^ 1 ', j^ 1 '), that is, contributions from each individual nucleon (for deuteron, 
there are two), include the charge density: 

p( 1 )=e^i±^(2^ 3 )(fc + ^- Pj ), (1) 

i=l 

and the 3-current densities from spin and motion of nucleons: 

2 i 

js P in = eV^« x ( p i % - Pt )(2Kf5^{k + p' t - Pi ), (2) 
2m N 

i=i 

jconv = eY,^^{p' l +Pi){^f^ ) {k + p' l -p l ), (3) 

i=l 

in momentum space, where fx l N is defined as 

with [i s = 0.88 and fi v = 4.71; and fc, p'^ and pi denote the 3-momentum of outgoing photon, 
outgoing ith nucleon, and incoming ith nucleon, respectively. 

As the canonical picture of nucleon-nucleon (NN) interaction is realized through the exchange 
of mesons, the nuclear EM currents should also contain these exchange effects, which are two-body 
in character. Since we are still lacking of a fundamental understanding about the AW interaction, 
the exact exchange currents (ECs) are also unknown. However, in order to reduce theoretical 
uncertainties in the studies of nuclear EM processes, it is important, given a chosen model for 
the NN interaction, to construct these two-body currents which are constrained by both current 
conservation and phenomenology. 
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For the PC NN interaction, many high-quality model potentials exist, and our choice for cal- 
culation in this work is the Argonne v\g potential (Avis) [111]. Although it gives good fits to the 
scattering data and deuteron properties, it is not straightforward to construct the corresponding 
ECs because the connection with the meson exchange picture is not clear for some parts in this 
potential. One traditional way to construct the ECs is implementing the NR minimal coupling 
(MC) to the potential, i.e., 

p_p_£(l + r*)A; H^H+ e -(l + T z )A\ 

then identifying the EM currents from the interaction Hamiltonian density, ej^A^. However, this 
procedure only constrains the longitudinal components, while giving no information about the 
transverse components which are conserved by themselves. Some uncertainty about these last 
terms comes from potentials involving quadratic velocity-dependent components as discussed in 

ilef. Q|. Moreover, the derivation of exchange currents for a model like Auis, employed in Ref. 

131 ] for instance, requires further elaboration. This model contains a Gaussian type component 
while the above prescription is usually applied to Yukawa potentials. Therefore, we leave this as 
an open question for future work and follow the treatment of Ref. Q] to examine: (1) to what 
degree current conservation is broken by the omission of PC ECs, and (2) how much the inclusion 
of PC ECs due to the one pion exchange, which gives the long-range part in A^ig, could restore the 
conservation. By this exercise, one can get some qualitative handle on this problem. The detail is 
be to discussed in Section llVl 

For the PNC AW interaction, our choice is the potential suggested by Desplanques, Donoghue, 
and Holstein (DDH) jlfl| . Since this potential, based on the one-boson exchange scheme involving 
7r, p, and uj mesons, has a close tie with the exchange picture, a more field-theoretical formalism, 
the so-called 5-matrix approach 



nding ECs. As 

we will show later, some transverse components arise naturally in this derivation. For clarity, we 
divide the following discussion into three parts: first, the model Lagrangian, consistent with the 
DDH scheme, is constructed; second, the PNC ECs are derived; and finally, we show how these 
exchange currents allow one to fulfill current conservation, given the DDH potential. 
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A. Model Lagrangian 



The total Lagrangian density we consider is expressed as 



C = £ + £pc + £pnc + £em + 5C , (4) 

where 5C contains all the terms not relevant for this discussion. The free Lagrangian densities for 
nucleon, N; pion, rho, and omega mesons, 7T, p M , and are 

C = N*(# - m N )N + i(<9 M 7r) • - \mW - If$ • + \m 2 pP , ■ - ^(fyp") • ( V) 

_1 F H F H/- + Im^^M _ ±(^)(<W0 , (5) 

where and i 7 }^ are the field tensors of the vector mesons, and we keep the gauge-fixing 
terms of vector mesons explicit for the moment. The PC and PNC meson-nucleon interaction 
Lagrangian densities are 



£-PC = ■ignNNN'y 5 T ■ tvN - gp NN N'(^fj, - i-^-^a pv q u )T ■ p^N 

-9 U nnN'(^ - i^a pu q u )^N , (6) 
2m N 

Cpnc = "^pV x n) z N + N'[h°T ■ p» + h\f» + ^(3rV " r ■ p»)] W N 

+N'(h° ul ^ + hir z ^) W N, (7) 

where q 11 is the 4-momentum carried by the outgoing boson; the strong couplings, g X NN, as well as 
the weak couplings, hx , are of DDH's definition (except their PNC pion coupling, f n , is renamed 
as h\ here); and the anomalous strong isoscalar and isovector magnetic moments of nucleon, Xs 
and Xv, are assumed to be the same as the EM values, -0.12 and 3.70, by vector meson dominance. 
Now we apply the covariant MC: 

to the above Lagrangian densities to obtain the EM interactions. For our purpose, only terms of 
first-order in e are included in Cem- Prom Co, we get 



4m 7 ' = -eN'm S \Q 2 )\ + FPiQ 2 ) y )7„ " ^(^W)^ + ^>r<FF& >(») 
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Figure 1: The vertex factor for vector-meson-photon coupling, where a, /?, and 7 are the Lorentz indices; i 
and j are isospin indices. 

4m 7) = -e(7rx<Vr)^, (9) 



(PP7) 



'EM " ~<P U >< - \<P» >< V)^ • (10) 

Note that in order to account for the nucleon structure, nucleon EM form factors, (superscript 
"S" for isoscalar and "V" for isovector; subscript "1" for Dirac and "2" for Pauli) have to be added. 
At Q 2 = -q 2 = 0, F[ s) (0) = Fl v) (0) = 1, F 2 (S) (0) = -0.12, and F^ v) (0) = 3.70. In principle one 
should also take into account the meson structures; however, they are still poorly constrained so 
we simply assume them elementary. 

Due to the momentum-dependent coupling of p meson to the nucleon anomalous magnetic 
moment, a Kroll-Ruderman type contact interaction reading as 

4m P7> = -e^^N'a^r x p»fNA» , (11) 
2m N 

also arises. This will lead to a seagull current which is important for current conservation, but was 
ignored in Ref. P|. 

It is worthwhile to point out that the EM interactions obtained above depend on the model 
Lagrangians we start out. For example, compared with the result of QHD II [l^J, a gauge theory 
approach which first models pion from a global SU (2) symmetry group and then adds p meson and 
photon by a local SU(2) ® U(l) symmetry group (both steps involve Higgs mechanism to generate 
meson masses), one observes a larger pfry interaction in QHD II by the amount of 

AZ;^ = \e(p, x p v fF^ , 

which modifies the pp-y vertex in an interesting way as we are going to explain. 

By fixing the gauge parameter £ = 1, i.e., 't Hooft-Feynman gauge, we obtained the same C^^f 



as in Ref. |lj|, and this gives a vertex factor (see Fig. QJ: 

esy l(qi ~ Q2Yg aP + k a g^ - fcV°] , (12) 
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with k + qi + q2 = 0. The latter two terms are combined to give amplitudes conserved by themselves 
when dotted by k^, which means they are purely transverse (actually, they correspond to magnetic 
dipole couplings which explain the transversality) . By adding AC^]\ it simply doubles the self- 
conserved terms, so that we have 

e 3ij [(<Zl ~ q2T9 aP + 2fc V - 2fc Vl , (13) 

for the vertex. As p meson is a spin-one particle, it has charge (c), magnetic dipole (ji), and charge 
quadrupole (Q) couplings to the EM field. When assuming it is an elementary particle, c = e, 
fj, = e/rrip, and Q = — e/m 2 , the vertex factor appears to be Eq. (fTTTf) 0,01) so the MC result 
under-predicts the p meson magnetic moment by two. This factor of two difference in self- conserved 
terms between MC and chiral Lagrangian approaches has been pointed out in Refs. 3, 13| It 
is also obtained from a quark model calculation of the charged p meson magnetic moment (in the 
limit where the p mass is taken as twice the constituent quark one). Although the factor two 
originates from different models, in order to have a closer contact with phenomenology, we make 
this particular modification to the ppj vertex. 

The modification mentioned above is just an example of model-dependences in constructing ECs. 
Since these purely transverse terms, often called non-Born (NB) terms, could not be constrained by 
current conservation, it is not easy to set up criteria a priori to judge which should be included or 
not, unless compared with experiments . For our calculation, we choose to include pirj and unrj 
interactions 

4 P M 7) = e f^^^ (7W (P 7 • , (14) 
= ef^e Q/375 F(^V3V) , (15) 

where the total anti-symmetric tensor is defined as £0123 = — 1 Q, Q]> because these coupling 
constants could be determined from the decay data (except for signs); and we ignore all the nucleon 
isobaric excitations for they are not the main theoretical emphasis of this work. The present work 
could be easily extended to them if necesssary. 

To sum up, the total EM Lagrangian density we consider is 

£ EM = 4m 7) + 4m 7) + {l { zT + \<p, >< p») z f^} + 4l7 ) + 4m 7) • (16) 



1 There are a lot of discussions in literature addressing this issue. We refer readers to more extensive reviews such 
as Refs. 0, (and references therein) for details. 
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Figure 2: Classification of meson exchange currents: (a) norm-recoil, (b) pair, (c) mesonic, (d) seagull, (e) 
isobaric, and (f) non-Born mesonic, where N and N* denote nucleon and nucleon excited state; M and M' 
denote mesons. 

B. PNC Meson Exchange Currents 



Diagrammatically, the ECs could be classified, according to Fig. as: (a) norm-recoil, (b) 
pair, (c) mesonic, (d) seagull, (e) isobaric, and (f) NB mesonic types. The division of norm-recoil 
and pair terms simply comes from the separation of positive- and negative-energy components in 
the covariant nucleon propagator. Sometimes confusion arises when it comes to pair and seagull 
diagram. If the PC irNN coupling is formulated as pseudo- vector, the seagull term is 0(l/m N ), 
while the pair term is higher-order in l/m N . On the other hand, if the pseudo-scalar coupling is 
adopted as we do here, there is no seagull term, however, at the leading order, the pair term looks 
exactly as the seagull term in the pseudo- vector scheme. Therefore, as far as the NR approximation 
is valid, these two formalisms are equivalent For the case of p meson, both pair and seagull 

diagrams have 0(l/m N ) contributions. As for the NB contributions from (e) and (f), we only 
consider the latter as explained above. 

Before one applies Feynman rules to evaluate these diagrams and extract the corresponding ECs, 
the gauge parameter has to be fixed. Though the physical results should be gauge-independent, a 
proper choice may greatly simplify the calculation. Here, we adopt the 't Hooft-Feynman gauge, 
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£ = 1, for the following reasons. First, the propagator is simpler 

<vu = 2 q ; qu 2M . ) 

\ q z — t,m z + te J 



q 2 — m 2 + it 



— m z + ze 

Second, the PNC potential, constructed from NR reduction of the one-boson exchange diagrams, 
corresponds to the form given by DDH. The last and most important of all, in combination with the 
previous statement, the contribution from the norm-recoil diagram represents how the one-body 
EM matrix element is modified by the presence of this NR potential Therefore, it should not 
be double-counted if one has already taken care of this by using the perturbed wave function, the 
route we will follow. 

In momentum space, to the order oil/m N , we have the following results for the pair and p-seagull 
(KR) 3-currents 

= ^7i^r (T1 ' " TiT2)[ai] 4 + (1 2) ' (17) 

Jpair+KR = [ (h° P (n • r 2 + r|) + -%(3rfr| - Tl • r 2 + [-1 - era] 



ipair 



+fcj(l + rf )[^ - rf<7 2 ] + (1 + Xv) K - (r, x r 2 )>, x «r 2 ] 1 

x P!W0 + (1 „ a) , (18) 

<?2 + m p 

-ew (1 + rf)(/l o [(7i _ ^ + [TiVi _ r|(72]) (2ttW_0 + (1 ^ 2) (19) 



2M/v <?2 + m u 

where the p-seagull 3-current corresponds to the term involving Xv', Qi,2 = v'\ 2 ~~ Pi, 2; and the 5 
function imposes the total 3-momentum conservation, k + q\ + qi = 0. All the pair charges are of 
higher order in l/m N compared with the nucleon charge, which is O(l), so they are neglected. 
For the mesonic 3-currents, we obtain 

-e9n NN K,_ .... , „ {2tt)HW{ 



Tl-T2- T[Tl)[q 2 - qi}a\ ■ giT-g- 2W 2 , 2T + (1^2), 



+2[(p' 1 + pi) + i(l + Xv)^i x qi]<T 2 • fc + 2[<x 2 ] {2m N k - k ■ ((p[ + p x ) + i(l + Xv)<r\ x 

(2vr) 3 ,5( 3 )( 



(<?i + m2 P )(<l2 + m l) 



+ (l<-2), 
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'mesonic 



imesomc 



e-gpNNQfm^hl (2vr) 3 5( 3 >(- • • ) 

-(n x t 2 ) 2 [qi x q 2 ] _ 2 \(„2 , 2\ + (* ^ 2 ) , 



(q( + m 2 )(q£ + ml 



0. 



where ko = E\ + E% — E[ — E' 2 = E} — Ef. Specially note that there is a contribution to the charge 
density at the same order as the nucleon charge: 

" "2e W ( K - i( Tl x r 2 )V 2 • k , W* 6 ^ ■ 



mesonic 



2^ 



( q ?+m2)(g2 +m 2 ) +( 1 - 2 )- (M) 



C. PNC NN Interaction and Current Conservation 

By considering the one-boson exchange diagrams where one of the meson-nucleon couplings is 
PC and the other PNC, the DDH potential in momentum space could be expressed as 



yn 

PNC 



v p 

' PNC 



where 



2\/2m N 

~9pNN 



i (n X T 2 ) 2 (<Ti + <T 2 ) • U n , 



h° p TL ■ t 2 + ^(7f + r|) + ^(Srfrl - r x ■ r 2 ) 



X (i(l + Xv) (<Tl X <T 2 ) • M p + (<T! - <T 2 ) • V p ) - - TI)(CT X + <T 2 ) • U 



QuJNN 

m N 



h l 

hi + -f(rf + r|) ) (i(l + Xs)(o-! x cr 2 ) • u, + (<n - <x 2 ) • w w ) 



+f(r!-r!)(a 1+ a 2 ).v u 



u 



(p[+Pi) - {P2+P2) 



x ~ 2 1 2 ' x 
q z + 



(25) 



(26) 



(27) 



2 (q 2 + m\) 

with q denotes the meson 3-momentum. 

To prove the conservation of these PNC currents at the operator level, we showed explicitly the 
following matrix element identities (with bra (Pi,p 2 | and ket |pi,p 2 )): 



([// \ V PNC \) — k • (jp a i r + j 'mesonic) ; (28) 

{[P^ j VpNc]) = " (jpair+KR + 3 mesonic) ' (29) 

([P (1) ,^vc]> = *•(#*■>. (30) 

mesonic ' mesonic ), (31) 
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u \Jmesonic 1 J mesomc/ ) V" / 

where if is the total Hamiltonian, which is the sum of kinetic energy (T) and both PC and PNC 
potentials (V PC and and the p mesonic current, Eq. I|21jh is separated into two parts: (I) is 

proportional to the vector (q\ — 92), and (II) contains the rest. The continuity equality of Eq. (|3*T]) 
indicates that (p meson i c , Jmellnic) forms a conserved current not constrained by the DDH potential, 

which actually results from the self-conserved p/ry vertex mentioned above, while the last equality 

(2) (2) 

shows the transversality of NB currents. Obviously, the total PNC EM EC operator, (p P M C > Jpnc) 
we construct, Eqs. I|17I24|) . satisfies the total current conservation condition 



[p« + pfl , T + V PC + V PNC ] = k • (j« + M + jfl 

(2) .(2) 

as long as the two-body PC EC, \Ppci3pc)-> 18 conserved, i.e., 



[p« , V PC ] =k j 



:(2) 



PC ) 



(2) 

(Ppc is higher order in l/m N ). Therefore, at least for the PNC part, we have every conservation 
condition met. 

For the calculation of AM, we have to use both the DDH potential and current operators in 
coordinate space, these expressions could be found in Appendix El 



III. DEUTERON ANAPOLE MOMENT 



Determination of the Deuteron Wave Function 



Due to the PNC AW interaction, the deuteron wave function, mainly a 3 S\ state with some 
fraction of 3 D\ component, could have parity admixtures in 3 Pi and 1 P\ channels. The former 
channel, induced by the isovector part of DDH potential, is dominated by 7r-exchange, while the 
latter one, resulting from the isoscalar interaction, is purely through the heavy-meson exchange. 
Therefore, we express the full parity-admixed deuteron wave function as 
1 



Mr) 



4nr 



<r) + -^«;(r) ) Coo 



' 3 \/3 
\l -{&! + <r 2 ) • rv 3p i(r) Cio + i-z- (fi - o" 2 ) • rv lp i(r) Coo 



(33) 



where S\2(r) = 3<ti • r <T2 • r — ct\ ■ <T2 and x an d C represents spinor and isospinor respectively. 
The Schrodinger equation in the center of mass frame is 



1 (\ 8 2 



1(1 + 1) 



m N \r dr 2 



V c {r) + V T {r) S 12 (f) + V PNC (r) 



*(r) 
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(34) 



where Vc(r) and Vr{r) are the central and tensor parts of the strong potential, respectively, and 
VpNc( r ) is the sum of the PNC potentials given above. Up to the linear order of PNC NN 
interaction, radial wave functions of each channel satisfy the differential equations 



u"{r) + m N (E — Vc{r)) u(r) = V8m N Vr(r)w(r) 



w 



'(r) - -yii)(r) -m N (E- V c {r) + 2V T {r))w{r) = VSm N V T {r)u{r) 



(35) 
(36) 



L(r) + ^=w(r)j ^ (i£(r) + V2F}(r) - V2F l J 



^SpiW " ^2^i(r) + m N (E - V c {r) - 2V T (r))d 3pl (r) 
2 

+2V2 (*J(r) - f*(r)) |- (u(r) + ^«i(r)) - ^ (*» - f*(r)) (u(r) - V^r 
«i p i - ^2^ipi( r ) + rn N (E - V c {r) + &V T (r))vi pl (r) = 



(37) 



2 

7! 



(u(r) - V2w(r)) {3 Xv F° p (r) - Xs F°(r)) 



-2 (3F p °(r) - F°(r)) — (u(r) - y/2w(r)) + - (3F°(r) - F°(r)) (u(r) + 2^w;(r) 



(38) 



where i^(r) = g„NNhlU(r), F p °(r) = g p NNh° p f p (r), F p l (r) = g p NNh l p f p {r), F°(r) = g U NNhlf u (r), 
and F^(r) = guNNh^f^r). In our numerical calculations, g nNN = 13.45, g pNN = 2.79, g U NN — 8.37, 
as well as DDH best values (in units of 1(T 7 ): h\ = 4.6, h° p = -11.4, h l p = -0.2, h 2 p = -9.5, 
/i° = —1.9, and = —1.1, are assumed. 

B. Anapole Moment : Expressions of Matrix Elements 



The anapole operator we use takes the form 

a = J dx x x (a? x jX 33 )) 1 (39) 
where j(a?) contains all the currents discussed in Sec. [HI Note that this form is equivalent to what 



lave been recommended in Refs. [9j,l25j, a result from implementing the extended Siegert's theorem 



TT 



With the deuteron wave function, Eq. (|33p . we obtain the anapole moment for the spin term; 

fj, v J drr (u(r) - V2w(r)j v 3p i(r) 
-V2~n s J drr (u(r) + -^=w(r)\v lp i(r) el, (40) 



a 



spin 



V6 
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where / = (°"l + a 2)xij z 1S the intrinsic spin taken in the spinor basis, and we note that 

this is equivalent to the total angular momentum taken in the total angular momentum basis, i.e., 
1 = (J = l, J Z \J\J= 1,J Z ) 2 . 

The matrix element of the convection current is written as 



3conv(x) — j^o nv {x) + j > conv (x) , 

Si 2 (r) 



ii S dridr2 7^rA j > r (r) + ^r w(r) ~ ^ (<T1 ~ a2) ' ™ lpl(r) 



x(Pi, P2V u(r) + 



a/3 \ 1 

^_ w{r)+i—{a 1 -a 2 )-rvi pl (r)\-j=- X iJ z , 



(41) 



Jconv ("") 



4m. 



Sl 2 (f 



u(r) + 



V8 



^(0 (Pi, P2) -*V o + cr 2) -ru 3 pi(r) 



-(o-i + <r 2 ) • ru 3 pi(r)(pi, p 2 ) I + w(r) 



4-7rr 



where 



(Pi, P2) ± = {pi, ^ 3 (« - ri)} ± {p 2 , ^ 3 (a; - r 2 )} • 



As been argued in 



181 



2vr 
~3~ 



7TC:' 



12m, 



(42) 



(43) 



(44) 



where Z = r x p 3 . Indeed, it can be shown that the matrix element of the operator (pi, p 2 ) + is 
proportional to 



R p+ -r P 



i? 2 + ^r 2 ) 



where R = (ri + r 2 )/2 is the coordinate of the center of mass and P = pi + p 2 is its conjugate 
momentum. In the center of mass frame where only the relative coordinate or momentum are 
relevant, the terms proportional to R or P are discarded. Consequently, the contribution from 
(Pi) P2) 4 " is zero and the magnitude from the convection current is determined by Eq. I|42jl. The 
anapole moment contributed by the convection term then reads 

1 7T 



a, 



3^6 



m» 



drr \u(r) — ■\/2w(r)j v^ p i(r)el . 



(45) 



2 Since the anapole moment is a vector moment, it should be proportional to the total angular momentum which is 
the only intrinsic vector of a system. 

3 We stress that the equality holds for the matrix elements and the operators themselves do not satisfy the equality. 
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Contributions from the PNC exchange currents are evaluated with the parity-even channels ( 3 Si 
and 3 Di) in the initial and final states. These channels are in spin-triplet, S = 1, and isospin-singlet, 
T = 0, state and for these states we have the spin and isospin selections rules 

{S = l\\{a 1 -a 2 )\\S = l) = (S = 1\\((t 1 X(t 2 )\\S = 1)=0, (46) 

{t = o\\(t 1 xt 2 ) z \\t = o) = o, (47) 

where the double bar 1 1 means the reduced matrix element. It is easily seen that the matrix elements 
of jpair+KR> jCesonid Jmesonici and the P art proportional to <t x - a 2 in j p pair and j£ air vanish. Then 
we are left with 7r-pair, 7r-mesonic, and the parts proportional to <T\ + a 2 in p°-pair and w-pair 
terms. Pair and mesonic pion contributions were already calculated in which are given as 



9m 



N 



2 



<air = ~ V Z'"" / drr'Ur) L(r) + ^=w(r) \ (u(r) - V2w(rj) e I , (48) 



V^Trg^NN f ( 1 

arr j n {r) I u(r) H — -=w{r) 



3m N m n J n \ 

u(r) ^1 - - -^ w ( r ) f 1 + ^ m 



and p°- and w-meson contributions read 



p° 27T g pNN 

~ pair ~~ 9m N 



<air = - Z 7^f J drr 2 f p (r) (u(r) + (u(r) - y/*w(r)) elh 1 ^ (50) 

2vr Qujjsin 



apair ' 9m, 



J dr r 2 f u (r) (u(r) + -^=w{r)^j (u(r) - y/2w(r)^ el hi,. (51) 



C. Numerical Results for the Argonne V\% NN Interaction Model 

With the Avis model, we obtain the numerical results (in units of fm 2 ) 

a spin = -0.547 /ii el+ (-3.7&J + 10.8 hi + 7.3 h° p + 3.7 h% ) x l(n 3 eJ , (52) 

a conv = 0.039 h\el+ (2.7 h\ - 7.6 hi) x 10~ 4 e I , (53) 

< = <ur + C« = ("0.027 + 0.028) h\el , (54) 

< = a£ ir = -0.7xlO- 4 fcJ e J, (55) 

< = a- mr = 1.8xl0- 4 /iie/. (56) 

The contribution from the nucleonic AM, a^r, should also be included in the full result. Since 
the deuteron is an isosinglet state, only the isoscalar component of the nucleonic AM contributes, 
i.e., 



2 

a N = (d\ ' 

i=l 



]T (o? } + a?> t?) ^ |d) = 2 a« (1 - |p D ) J , (57) 
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where a^y denote the isoscalar and isovector nucleonic AM, and P D is the probability of deuteron 



, y, y, y, q, and 



D-state. Several theoretical estimates for the nucleonic AM exist [2£ 

here we use the result of Ref. [^jj because it is the most recent one which includes the full E 
interaction at the nucleon level. 

a 



For the pion sector, Ref. |3JJ] gave 



7T(1) 



9 A h\ „ _ 6 m n ^ A x 

7r m N m n 



where A x is the chiral symmetry breaking scale (the authors chose it to be 47r/,r) 4 . When setting 



A x = m N , the leading term is equivalent to what has been used in Refs 



for the deuteron 



AM calculations, while the full result is the same as been used in Ref. pj. By including the heavy 
mesons, the numerical result is |3l| 



a (1) 



-0.274 /i* - 0.419 hi - 0.129 h° , 



(58) 



and this leads to the nucleonic contribution 



a N = (-0.250 hi - 0.383 hi - 0.118 hi ) e I . 



(59) 



Although we do not have the exact PC EC corresponding to Av±s potential, we try to ap- 
proximate this by the EC coming out from one-pion exchange diagrams. Corresponding currents 
in the configuration space are given in the Appendix A. Pair- and mesonic-term of the PC EC 
contributions read 

2 



a 



PC 
pair 



a 



PC 

mesonic 



1 ( th 
1 fth 



dre m ^ r (1 + m^r) v 3p \{r) (^u(r) — ^w^^j el, 



3^ V. 2m N 
Numerically we obtain 



dre m * r v 3p i(r) 



3u(r) — m n r ( u(r) + — ^=w(r) 



(60) 



el. (61) 



a 



PC 



a p c +a Pc 

"'pair ' ""mesontc 



-(7.8 hl + 0A hi - 1.1 hi) x 10~ 4 el. 



(62) 



In the next section, we will discuss if this is a reasonable approximation and to what extent is the 
breaking of current conservation due to the non-exact PC EC. 



4 Note that a factor of <in/m 2 N is needed for converting the dimensionless AM defined in Lorentz-Heaviside units 
into our definition. 
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Figure 3: The left panel shows the behavior of r 2 (u(r) + w(r)/^/~2)(u(r) — \/2uj(r)). The central panel 
compares the Yukawa function of the pion (solid line) and rho- meson (dotted line) . The right panel is the 
behavior of the integrand for the it— and p— pair terms in Eq. i!48|) and Eq. J 5011 , respectively 

Finally, the full deuteron AM could be expressed as: 

PNC PC 
dd — O'spin (Lconv ~i~ ^ex ^ex 

= (-0.756 hi- 0.387 h p + 0.010 hi + 0.007 h°- 0.114 h°jel. (63) 
IV. DISCUSSIONS 

Presenting the results, we separated the contributions of the pion and the heavy mesons so that 
their relative order of magnitude can be compared easily. In the result of the nuclear part, i.e. the 
contributions from deuteron wave function or exchange currents, which include spin, convection, 
PNC EC and PC EC, the magnitudes of the heavy meson contributions are commonly suppressed 
by an order of 2 or 3. This suppression can be understood most easily from Eqs. I|48 [ I50 [ IBT jl . PNC 
pair terms of ir, p and u mesons, respectively. Aside from the weak coupling constants, they differ 
only by the Yukawa functions f^{r) or f P)UJ (r) in the integrand. At small r, the deuteron wave 
function is proportional to r (for the simplicity of argument, we neglect D— state wave function) 
and increases very fast up to r ~ 2 fm at which the maximum is reached. Afterwards, the wave 
function decreases very slowly, converging to zero. The factor r 2 multiplied by u 2 (r) gives rise to 
the suppression at short range. The left panel of Fig. 01 shows this behavior. One can expect that 
this r 4 behavior at short distances makes the contribution from r < 1 fm quite small. On the 
contrary, the Yukawa function, which behaves like 1/r at small r and also depends strongly on the 
mass of the meson, is short-range peaked. A comparison of /^(r) and f p (r) is given in the central 
panel of Fig. [3 The quantity f p (r) is non-negligible for about r < 0.5 fm but the remaining part 
of the integrand, approximately r 2 u 2 (r) is small in this region. Consequently the heavy meson 
will be suppressed substantially compared to n. The right panel of Fig. 01 shows the behavior 
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term 


spin 


conv. 


PNC EC 


PC EC 


nucleonic 


total 


(p + U))/lT 


4.0 


0.6 


-2.0 


-3.1 


-26.2 


-5.9 



Table I: Comparison of the heavy meson with the pion contribution to the deuteron anapole moment, term 
by term. The ratios are given in %. 



of the total integrands in Eqs. If we approximate the integrand as r 3 exp(— m x r), the 

maximum of the integrand is reached at r ~ 3/m x . This approximation is good for p but bad for 7T, 
however, this can provide a semi-quantitative explanation for the suppression of the heavy mesons. 
At the maximum, the integrand has the value proportional to l/m x . If comparison is made for 7r 
and p, the maximum value of p is smaller than n by a factor of 1/200, which can account for the 
suppression of order 2 of the heavy mesons. 

To a lesser extent, a similar argument can be applied to the qualitative understanding of the 
heavy-meson suppression of spin and convection terms. The right hand side of Eqs. (|371 l38|) 
which are the sources of parity admixed P— states also contains the Yukawa functions. When the 
equations are solved, i.e. integrated with respect to r, one can expect some amount of suppression 
for the heavy mesons as was discussed in the analysis of pair EC's, but without the factor r 2 . In 
Tab. HI we summarize the ratio of the heavy meson contribution to the pion with the weak coupling 
constants given by DDH best values. The magnitudes of the heavy meson terms are suppressed by 
an order of two compared to the pion terms in common. As far as the nuclear part is concerned, tt 
contribution is by far dominant. If one can disentangle the nuclear part contribution from the total 
deuteron anapole moment, this can provide information to extract the magnitude of h\ with high 
accuracy. As a side remark, it is noticed that the combination of the isovector couplings, /i* , 
and h]j, entering the spin and convection contributions is close to the one determining low energy 
PNC nucleon-nucleon scattering, roue hly 3/4 + 0.02 /ii-0.06 hi Q. This lets us think that the 
Danilov's approach jjjjj, whose application for the deuteron anapole moment was proposed later on 
by Savage [6J, should provide a reasonable estimate for this contribution. 

Contrary to the nuclear part, the nucleonic one has sizable contribution from the heavy mesons. 
The relative ratio of the heavy meson is about 26% of the pion based on the nucleonic anapole 
estimate of Ref. j^l]]. However, we should also note that these authors considered more non-DDH 
type couplings such as non- Yukawa type irNN couplings and the inclusion of hyperons. So far, 
no detailed knowledge about these exotic couplings exist, and the theoretical uncertainty could be 
huge. The study of these terms and their implication for two-body nuclear contributions will be an 
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interesting topic for further exploration. 

Another issue that should be addressed is the gauge invariance of the results. We showed in 
Sect. [H]that the PNC ECs we constructed satisfy the gauge invariance with the DDH potential. 
However, with the phenomenological strong interaction models like the one adopted in this work, 
gauge invariance may not be satisfied if the ECs are not consistent with the phenomenological 
potentials. For Av\%, the potential contains 18 types of operators but we took into account the 
dominant ECs only, which naturally leads to inconsistency, or breakdown of gauge invariance. Even 
though the inconsistency may only change the result slightly, the investigation of the extent to which 
gauge invariance is broken could help to get an estimate of the error. 

Most phenomenological potentials are very complicated and the analytic analysis of gauge invari- 
ance is a formidable or practically impossible work. However as suggested in [8J, one can estimate 
the amount of gauge invariance breaking by comparing the anapole moment obtained from Eq. (j3*9"|) 
and 



a = —ir J dxx 2 j(x) , 
which can be obtained from Eq. with the assumption 



(64) 



V- j(x) = 0. 



(65) 



1*1 



In [a], we observed that the inclusion of the PC ECs of tt removes almost all the inconsistency. In 
this work, we follow the same procedure to investigate the breakdown of gauge invariance. 

The spin current can be shown easily to satisfy Eq. ijfiHjl . The convection, pair and mesonic 
anapole moment of tt, p— and to— mesons with the definition given by Eq. (|64*)l can be calculated 
straightforwardly to give 



TT 



a, 



pair 

1T,CC 

^mesonic 



a. 



2\^6m N 

ft 9ir NN 

2y2m N 

9-knn 



dr r z 



u(r) [ v 3pl (r) + 2 — - — ) + [ v 3pl (r) y — -L 



V2 



24^/2 



drr f n {r) (r) 



w 2 (r) 



hi el, 



e 7,(66) 
(67) 



u 2 (r) 



w 2 (r) 



1 



a 



a 



p°,CC 
pair 

u>,CC 



, \ w ( r ) 

(jn^r + 4) _ _ ^ u (r) + -A^ 



el, 



m n r (m n r + 3) 



pair 



r irr 2 Ur) ( u 2 {r) _™M) ti:eI 



2 m> 



hi el, (68) 

(69) 
(70) 
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Eq. !p) 


Eq. (|Pt 




K 


hi 


K 


K 


hi 


conv. 


0.03925 


0.00027 


-0.00076 


0.05348 


0.00024 


-0.00066 


PNC 7r pair 


-0.02668 






-0.07895 






PNC 7r mesonic 


0.02830 






0.04706 






PNC p pair 




-0.00007 






-0.00023 




PNC oj pair 






0.00018 






0.00061 


PC ir pair 


0.01400 


0.00013 


-0.00038 


0.06325 


0.00062 


-0.00177 


PC 7r mesonic 


-0.01478 


-0.00017 


0.00049 


-0.04366 


-0.00047 


0.00133 


Total-PC 


0.04087 


0.00020 


-0.00058 


0.02159 


0.00001 


-0.00005 


Total 


0.04010 


0.00016 


-0.00047 


0.04117 


0.00017 


-0.00049 



Table II: Coefficients of the weak coupling constants for given terms with different definitions of the anapole 
operator, Eq. (j^HJ and Eq. l(64*j) . 



= ^ (f^) / dre-^ r (l + m7T r)v 3pl (r) (u(r) - u,(r)) el, (71) 
12^/6 \2'> 



a mesonic = ~T7TTf WZT~ / Ar e u 3pl( r ) 



1 I 9-KNN \ I ^ rr - 

n N J J 

u(r)(18 + Im^r — m^r 2 ) jj- m n r (4 + m n r 



w(r) 



el, (72) 



where the superscript CC denotes the quantity calculated with Eq. (|64"|) . Numerical results are 
summarized in Tab. [D] 

Judging from the "Total" column in Table II, two definitions of the anapole operator only 
differ in results by 3%, 6%, and 4% for the 7T-, p-, and w-exchange respectively; this means our 
result satisfy current conservation very well. This is quite surprising because only the PC one-pion 
exchange current — not fully consistent with the adopted strong potential, Av±s — is included in our 
calculation. The reason is partly because the AM is a r 2 -weighted moment, plus the deuteron wave 
function peaks around 2 fm with a long tail; therefore the long-range physics, which is dominated 
by the one-pion exchange including the case of Avis, becomes much more important. This explains 
why our conclusion differs from the one found in Ref. where large effects due to undetermined 
ECs were found. Their observation was concerning a different operator (El transition), which has 
a range relatively shorter than the quantity considered in this work. 

Furthermore, one can observe that for the contributions from exchange currents, Eq. l(39|) always 
give smaller values than Eq. l(f)4"|) . This implies the calculation using Eq. (|39|l suffers less from 
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the incomplete knowledge or uncertainty of exchange currents which causes the breaking of current 
conservation. For example, if the PC pion ECs are left out in our calculation, by comparing the 
"Total" and "Total-PC" columns in Table II, the error is 2%, 25%, and 23% for vr, p, and u 
component respectively when Eq. (|3ll|) is used; while the error is -46%, -94%, and -90% for Eq. 

Perhaps, the reason for the small contribution of ECs in the former case is to be looked for 
in the proportionality of the dominant ECs to the position vector x, which readily gives zero when 
inserted into the corresponding anapole operator, Eq. (jHEJ- Therefore the merit of Eq. lj3*9"|) is 
justified, at least for the case of deuteron. 

In conclusion, we constructed the PNC ECs due to one tt— , p— and to— exchange, and showed 
that they satisfy the current conservation; and hence are fully consistent with the adopted DDH 
PNC potential. An application was made to the calculation of the deuteron anapole moment. We 
observed that, for the nuclear part, the contribution of heavy mesons is suppressed by an order of two 
compared to the pion, a result consistent with the similar work by Blunden Q|. Consistency with 
Av±8 was also checked. We found that the approximation of using only the PC one-pion exchange 
current is pretty good — the breaking of current conservation only amounts to a few percent, and 
this should be fixed by using the PC exchange current fully consistent with Avi&. Therefore, the 
contribution from the nuclear part to the deuteron anapole moment can be determined with an 
error less than 5%, while the major uncertainty should come from the nucleonic anapole moment 
instead. 
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Appendix A: PNC NN POTENTIAL AND EXCHANGE CURRENTS IN COORDINATE 
SPACE 

The DDH potential in coordinate space could easily be obtained by applying the following 
transformation rules to Eqs. 1)2 5127}! 

u x u x (r) = [p, f x (r)} , 
v x -> v x {r) = {p, f x (r)}, 
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where r = r\ — r 2 ; r = \r\; p = (p\ — p 2 )/2 = — «V r ; and the Yukawa functions f x (r) are defined 
as 

n — nix r 

fxir) 



4trr 

For the PNC MECs, we list all the leading-order, 0(1/777^), 3-currents, which are relevant for 
the AM calculation: 



j; air (x; r 1; r 2 ) = - e !h ™ K (n • r 2 - Tf r|) U(r) £ ^(z - r*) cr< 

x; ri, r 2 ) = -—7= (n • r 2 - r{ r|) 

2 y2 ni N 



(Al) 



imesonic 



x (Vi - V 2 ) [(o-i ■ Vi - tr 2 ■ V 2 ) , U(r xl ) U(r x2 )} , (A2) 

x (Vi • r 2 - Tfrf ) - tr 2 ) (tf( 3 )(x - n) - ^(sn - r 2 )) 

+(l + Xy)(TiXT 2 ) z (o- 1 xo- 2 ) ^(^(as-r^ , (A3) 



^ + g ^ C 7 ? + T 2 ) + ^= (*x ~ <*i) + - h p (rf - r|) (<n + <x 2 )J 

x ((i + 7 f),5(3)(a;-r 1 )-(l + 7f)<y®(a!-r2)), (A4) 
JSLntote ^ r 2 ) = - ^ ^ fa x r 2 y 



x|(Vx- V 2 ) ^-^^-^-{(Vi-Va), f P (r xl ) f p (r x2 )} 

+ (1 + X v) (o-i x a 2 ) ■ [(Vi - V 2 ) , f p (r xl ) f p {r x2 )\ \ 

+2 V a x (i {V? <x 2 - V a 2 <tx + a? V 2 - of Vi , / P (r x i) / p (r x2 )} 
-(1 + X v) [(ai x Vi) a er 2 - (a 2 x V 2 ) a o-i 

+of<7 2 x V 2 -ff2(TiX Vi, / P (r x i) / P (rx2)]l 

+4 7 777^ (c^ + <X 2 ) [if , / p (ra,i) / p (ra. 2 )] > (A5) 



Jpair-W n, r 2 ) = — — — Mr) 

Z 771 ?v 
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x f ( hi + \hi (rf + rf )) (<n - <r 2 ) + 1 (rf - rf ) (<n + <r 2 ) 

x ((1 + 7f ) <y< 3 > (a- - n) - (1 + r|) ^ (x - r 2 )) , 

(x; n, r 2 ) = _ e ^^7 ^ (ri x T2 y 
\/2m p 

x(Vi X V 2 ) (/ p (r s i) /tt^) + /tt(^i) /pOx2)), 



) 



(A6) 



mesomc 



(A7) 



where r^j = |x — Vj and V x act on the source point and field point x respectively; and the 
superscript a is the index to be summed from 1 to 3. Note that we separate charged and neutral 
p mesons according to their isospin structure; and the last term of Eq. (A5) should be combined 
with a charge density corresponding to Eq. (24) in order to insure the current conservation. 
PC one-pion ECs of the pair and mesonic terms that contribute to the AM read 



J^esonUx; n, r 2 ) = - (tl x r^Vi • Vi) (<r 2 ■ V 2 ) (Vi - V 2 ) U(r xl ) U(r x2 lA9) 
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